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Abstract 

The purpose of this note is to analyze the long term stability of the Navier- 
Stokes equations supplemented with the Coriolis force and the stress-free 
boundary condition. It is shown that, if the flow domain is axisymmetric, 
spurious stability behaviors can occur depending whether the Coriolis force 
is active or not. 



1. Introduction 

The liquid core of the Earth is often modeled as a heated conducting 
fluid enclosed between the solid inner core and the mantle. Numerically 
simulating the dynamics of the liquid core is difficult in many respects; one 
of the difficulties comes from the presence of viscous layers that develop at the 
boundaries of the fluid domain, i.e., the so-called inner core boundary (ICB) 
and core mantle boundary (CMB). It is a common practice in the geophysics 
literature to use stress-free boundary conditions in order to minimize the role 
played by the viscous layers. Although this choice of boundary condition is 
convenient, it is not clear that it is more physically justified than using the 
no-slip condition. Actually, enforcing either the no-slip or the stress-free 
boundary condition may lead to significantly different results when it comes 
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to simulating the geodynamo. For example, Glatzmaier and Roberts [5] and 
Kuang and Bloxham [8, 9] have used the above two different sets of boundary 
conditions and have reported numerical buoyancy-driven dynamos in rapidly 
rotating spherical shells that differ in some fundamental aspects, see e.g. 
[14]. The simulations reported in [8] use the stress-free condition whereas 
those reported in [5] use the no-slip condition. The dynamo simulated in 
[8] is composed of an external magnetic field dominated by an axial dipole 
component, like that of the Earth, with an intensity close to the present 
geomagnetic dipole moment. The external magnetic field is comparable to 
that obtained by Glatzmaier and Roberts [5], but important differences in 
the velocity and magnetic fields between these two dynamos can be observed 
within the outer core and the Taylor-Proudman tangent cylinder. (It is 
known that rotation of the Earth rigidities the flow field in the direction 
parallel to the rotation axis through a mechanism known as the Taylor- 
Proudman effect. This effect makes the imaginary cylinder that is tangent to 
the equator of the solid inner core and whose axis is parallel to the rotation 
axis of the Earth act like a solid boundary.) In the dynamo reported in 
[8] the fluid flow is almost stagnant inside the tangent cylinder and has a 
strong azimuthal component outside; the magnetic field is composed of two 
opposite toroidal cells and a simple dipolar poloidal structure and is active 
throughout the outer core. In the dynamo reported in [5] the fluid flow is 
composed of an intense polar vortex that is located inside the tangent cylinder 
and extends in the two hemispheres; the toroidal component of the magnetic 
field is active only inside the tangent cylinder and is concentrated near the 
ICB; the poloidal component has a complicated dipolar structure with extra- 
closed loops near the ICB. It is suggested in [14] that the significantly different 
structures of the above two dynamos should be attributed to the nature of 
the boundary conditions that are imposed at the ICB and CMB interfaces. 

In addition to thermal or compositional convection due to buoyancy, pre- 
cession is also believed to be a possible source of energy for the geodynamo. 
The precession hypothesis has been formulated for the first time in [1] and 
experimentally investigated using a water model in [11]. It has since then 
been actively investigated from the theoretical, experimental and numeri- 
cal perspectives. However, it seems that it is only recently that numerical 
examples of precession dynamos have been reported in spheres [15, 16], in 
spheroidal cavities [17] and in cylinders [13]. Recently, Wu and Roberts [17] 
have numerically studied the dynamo effect in a precessing oblate spheroid. 
To facilitate their analysis the authors have split the total velocity field into a 
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basic stationary analytic (polynomial) solution (the so-called Poincare flow) 
and a fluctuating part. Following ideas of Kerswell and Mason [12], they have 
implemented the stress-free boundary condition on the fluctuating compo- 
nent of the velocity in order to reduce the impact of the viscous layers at the 
rigid boundaries. 

The purpose of the present paper is to show that the use of the stress-free 
boundary condition poses mathematical difficulties. We prove for instance 
that, if the fluid domain is not axisymmetric, the flow always returns to rest 
for large times when the stress-free boundary condition is enforced, but this 
may not be the case when the flow domain is axisymmetric. Various scenarios 
can occur depending whether the domain undergoes precession or not. 

The note is organized as follows. We analyze the stress-free boundary con- 
dition in general fluid domains in §2. We show that this boundary condition 
is admissible if and only if the domain is not axisymmetric (see Proposi- 
tion 2.2). We revisit the same question in axisymmetric domains that un- 
dergo precession in §3 and §4. We show in §3 that the problem exhibits a 
spurious stability behavior if the stress-free condition is enforced on the ve- 
locity field minus the Poincare flow (i.e., on the perturbation to the Poincare 
flow). We show in §4 that the problem always returns to rest for large times 
if the homogeneous stress-free boundary condition is enforced. The theoret- 
ical argumentation developed in §3 and §4 is numerically illustrated in §5. 
Concluding remarks are reported in §6. 

2. Stress-free boundary condition without precession 

The objective of this section is to investigate the long term stability of the 
Navier-Stokes equations equipped with the stress-free boundary condition. 
The fluid domain is denoted Q and is assumed to be open, bounded and 
Lipschitz. 

2.1. Position of the problem 

We are interested in the motion of an incompressible fluid in a container 
Q with boundary T. The container is assumed to be at rest in a Galilean 
frame of reference. Denoting u the velocity of the fluid and p the pressure, 
the fluid motion is modeled by means of the incompressible Navier-Stokes 
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equations: 



d t u + u Vu - 2i/V-e(u) + Vp = 0, (2.1) 

Vu = 0, (2.2) 
U| t=0 = u , (2.3) 

where v is the kinematic viscosity, e(u) := | (Vu + Vu T ) is the strain rate 
tensor, and u is an initial data in H := {v G L 2 (f2) : Vv = 0, vn r = 0}. It 
is a common practice to replace the term V (Vu + Vu T ) in the momentum 
equation by Au since VVu T = for incompressible flows. We nevertheless 
keep the original form of the viscous stress since we want to enforce the 
so-called stress-free boundary condition: 

(n-e(u)) xn, r = 0, (2.4) 

together with the slip boundary condition: 

n-U| r = 0, (2.5) 

where n is the unit outward normal on T. The stress-free condition means 
that the tangent component of the stress at the boundary is zero. We shall see 
that this boundary condition is admissible in general for non-axisymmetric 
domains, but it yields pathological stability behaviors if the fluid domain is 
a solid of revolution. 

We are not going to discuss the well-posedness of the above problem in its 
full generality since it is still unknown whether the three-dimensional Navier- 
Stokes equations are well-posed under the much simpler no-slip boundary 
condition. We nevertheless recognize as a symptom of pathological stability 
behavior the fact that there are solutions to (2.1)-(2.2)-(2.3)-(2.4)-(2.5) that 
do not return to rest as t — > +oo if Q is axisymmetric. 

Definition 2.1. We say that Q is stress-free admissible if there is a constant 
K > 0, possibly depending on Q, so that the following holds 

K [ v 2 <: [ e(v):e(v), VvGH 1 ^), vn| r = (2.6) 
Jn Jn 

where ":" denotes the tensor double product. 

Proposition 2.1. Assume that Q is stress-free admissible, then {0} is the 
global attractor of (2.1)-(2.2)-(2.3)-(2.4)-(2.5). 
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Proof. We omit the details concerning the existence of Leray-Hopf solutions, 
which can be constructed using standard Galerkin techniques [10], and we 
focus only on the aspects of the question which are relevant to our discussion. 
It is clear that is an invariant set of (2.1)-(2.2)-(2.3)-(2.4)-(2.5). Let B be 
a bounded set in H and let u e B. Let u be a Leray-Hopf solution corre- 
sponding to the initial data u and let v be a smooth solenoidal vector field 
satisfying the slip boundary condition. Upon multiplying the momentum 
equation by v and integrating over the domain we obtain 

/ d t u-v+ / u-Vuv-2z/ / Ve(u) v+ / Vp-v = 0. 
Jn Jn Jn Jn 

Solenoidality and the slip boundary condition imply that f Q Vp-v = — f n pV-v+ 
/ r pv-n = 0. Now, using the decomposition 

v = (n-v)n — nx (nxv) , 

and integrating by parts the viscous term we obtain: 

- / Ve(u) v= / e(u):Vv- / n e(u) v 
Jn Jn Jt 

= / e(u):e(v)+ / (n-e(u)xn) -(nxv) = / e(u):e(v). 
Jn Jv Jn 

The transport term and the time derivative are re-written in the following 
form 

/ u-Vu v = / Jv(u(u-v)) + \ I (u-Vu v - u-Vv u) = J J (u-Vu v - u-Vv u) , 

Jn Jn% 2 J n 2 J n 

d t u-\ = \ ( <9 t (u-v) + \ I (d t u-v - d t v-u) . 



n % Jn 2 Jn 



We now apply the above identities by replacing v by a sequence {v n } neN 
that converges in the appropriate norm to u. By passing to the limit (we 
omit the details again), we finally obtain 



1 d 

2 ~dtj n 



u 2 + 2z/ / e(u):e(u) < 0. 
Jn 



Note that equality is lost in the passage to the limit. Whether equality holds 
in general is an open problem which is part of the Millenium prize. Then 
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using (2.6), we infer the following inequality: 



1 d 
2di 



J 

Jn 




which immediately leads to 



u llL 2 (n) ^ ||uo||L 2 (n)e 



-2Kvt 



thereby proving that u — > as t — > +oo. 



□ 



We shall see that the stress-free admissibility condition (2.6) does not 
hold for axisymmetric fluid domains, which are common in geoscience. 

2.2. The non- axisymmetric case 

To better understand the stress-free admissibility condition (2.6), we first 
prove that it holds if and only if Q is not axisymmetric. 

Definition 2.2. We say that fl is axisymmetric (or is solid of revolution) if 
and only if there is a rotation R : Q — > Q which is tangent on T. 

Upon introducing the average operator over Q, (v) : = p- f n v, where \Q\ 
is the volume of Q, the following lemma gives a characterization of non- 
axisymmetric domains: 

Lemma 2.1 (Desvillettes-Villani [3]). Assume that the domain Q is not a 
solid of revolution of class C 1 , then there is c> so that 



We are now in measure to state the main result of this section: 

Proposition 2.2. Assume that the domain fl is of class C 1 , then Q is stress- 
free admissible if and only if Q is not a solid of revolution. 

Proof. Let us assume first that Q is not a solid of revolution and (2.6) does 
not hold. We start from the Korn inequality (cf. e.g. [4]): there exists a 
constant c > such that, for all v e H 1 (fi), 



c|fi|(Vxv> 2 < ||e(v)|| L2(n) , Vv G H 1 ^), vn| r = 0. 



v|| L 2(fi) + ||Vv|| L 2 (n) < c (||v|| L 2 (n) + ||Vv + Vv T || L2(n) ) . 



(2.7) 



6 



Since (2.6) does not hold, for any nGN, one can find u„ G H 1 (fi) such that 

Tu 1 

u n -n| r = 0, ||u n || L 2 ( n) = 1, and || Vu n + Vu n || L 2 (n) ^ -. 

Th 

The Korn inequality implies that the sequence u n is bounded in H 1 (f2). 
Since the inclusion H 1 (r2) C L 2 (f2) is compact, there exists u e H 1 (f2) 
such that (we keep using u n after extraction of the converging sub-sequence) 
Il u n — u IIl 2 (si) - > and u n — ^ u in H 1 (fi). We also have 

Vu n + Vu^ ->■ in L 2 (Q) and Vu n + Vu^ ->■ Vu + Vu T in &(&,), 

which finally gives Vu+Vu T = (T>(Q) is the space of smooth vector- valued 
functions with compact support in Q and T>'(fl) is the space of vector- valued 
distributions over Q, i.e., the linear forms acting on T>(Q).) Applying the 
Korn inequality to u — u„ and using the fact that 

||u„ - u|| L 2 (n) + ||Vu„ + Vu^ - Vu - Vu T || L 2 (n) ->■ 0, 

we infer that ||u n — u|| H i(n) — > 0. This allows us to pass to the limit on the 
boundary condition u-ri|r = 0. The condition e(u) = implies that there 
are two vectors t e M 3 , w G I 3 so that u = t + wxx. This means that 
Vxu = (Vxu) — uj. Using Lemma (2.1), we conclude that uj — 0, which 
means that u = t. The boundary condition u-ri| r = implies t = 0; this in 
turn means u = 0, which is impossible because ||u|| L 2(q) = 1. In conclusion, 
(2.6) holds. 

Let us assume now that Q is axisymmetric. This means that there is a 
rotation R : Q — > Q which is tangent on T. Let us assume that the rotation 
axis is parallel to e z and the coordinate origin is located on this axis. Then 
R(x) = we 2 xx and clearly R e H 1 (fi), R(x)-n(x)| r = 0, ||R||L 2 (n) 7^ but 
(2.6) does not hold since e(R) = 0. □ 

2.3. The Axisymmetry curse 

Let us assume that Q is axisymmetric. We are going to show the following 
statement in this section. 

Claim 2.1. The zero velocity field, 0, is in the global attractor of (2.1)-(2.2)- 
(2.3)-(2.4)-(2.5) 7 but the rest state, {0} 7 is not an attractor. There are initial 
data that create flows that never return to rest. In particular, if the initial 
data is a solid rotation, the flow will rotate for ever without losing energy. 
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Recall that it can be shown that Q is axisymmetric if and only if Q is 
either a sphere (and all the directions are symmetry axes) or Q has a unique 
symmetry axis. Without a loss of generality, we assume Oz is the only 
symmetry axis of Q. Recall that all the solid rotations about Oz can be 
written as follows x i — > ue 2 xx, u G K, where x is the position vector. We 
introduce the following space 

K : = span{e z xx} (2.8) 

and its orthogonal in L 2 (f2), say 1Z ± . 

Lemma 2.2. Let Q be an open, bounded, connected, domain of class C 1 
with unique symmetry axis Oz. There exists K > such that, for every 
vG^n H 1 ^) with vn = 




where we denote |e(v)| 2 := e(v):e(v). 

Proof. The proof is similar to that of Proposition 2.2. By contradiction, we 
consider a sequence v n G TZ ± r\H 1 (Q) with vanishing normal component such 
that 

||v n || L 2 = 1 and ||e(v n ) T || L 2 ^ -. 

n 

Using Korn inequality, we can prove that (up to extraction) v n converges in 
H 1 (fi), and the limit v satisfies 

v G TZ L , v-n = and e(v) = 0. 

This implies that v is the sum of a translation plus a solid rotation. But 
Q being bounded the translation is zero and v is a solid rotation about 
Oz-axis (recall that f3 ^ 0), i.e., v G 1Z fl TZ L = {0}, which contradicts 
||v|| L 2(n) = 1- □ 

We claim that the Navier-Stokes problem (2.1)-(2.2) equipped with bound- 
ary conditions (2.4)-(2.5) has spurious stability properties due to the follow- 
ing proposition. 

Proposition 2.3. (i) K is the global attractor of (2.1)-(2.2)-(2.3)-(2.4)- 
(2.5). (ii) No element in 1Z is an attractor. 
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Proof, (i) Let u G L 2 ((0, +00); L 2 (ft)) n L°°((0, +00); H 1 (Q)) be a Leray- 
Hopf solution of (2.1)-(2.5) and consider the following decomposition: 

u(t) = u x (t) + A(t)e 2 xx, where u x (t) G 7^, A(t) el, Vt G [0, +00). 

u being a Leray-Hopf solution implies that 

l|u ± (t)||L(n ) + A(t) 2 ||e,xx||2 2(n) +4z/ / / |e(u)| 2 < ||u || 2 2(n) . 

Jo Jn 

which owing to Lemma 2.2 implies 

11^)11^) + ^)1e 2 xx|| 2 2(n) +4z/K f Hu^l^dr < ||u || 2 2(n) . 

Jo 

Using the Gronwall-Bellmann inequality, we infer that ||u J -(t) \\l 2 (q) < || u o||L 2 (n)e" 

Invoking Lemma 3.1 we infer that = 0, implying that A(f) = A(0). In 
conclusion 



-2vKt 



|u(t) - A e 2 xx|| L 2 (Q) = Hu- 1 ^)!!^^) < ||u ||L2(n)e" 



-2vKt 



This implies that the global attractor, say A, is such that A C 7Z, but since 
Ao spans R, we conclude that A = 71. 

(ii) Let us consider the solid rotation field u = ue 2 xx G 7Z. It is clear 
that u is invariant, i.e., is a steady-state solution. Let B(u, p) G H be 
the ball centered at u of arbitrary radius p > 0. Let v = pe z xx. G 7Z, 
p 7^ 0, be another solid rotation and assume that p is small enough so that 
u + v G B(u, p). Let us observe that 

(u + v)-V(u + v) = 2(u + v)-e(u + v)-(u + v)-(V(u + v)) T = -^V|u + v| 2 , 

since u + v is a solid rotation; moreover, u + v satisfies (2.2), (2.5) and 
e(u + v) = 0. The property e(u) = implies (2.4) and V(e(u)) = 0. Upon 
setting p=||u + v| 2 we conclude that u + v solves (2.1). This proves that 
u + v is invariant (i.e., a steady-state solution). In other words u + v does 
not converge to u, no matter how small p is, thereby proving that the set 
{u} is not an attractor, no matter how large v is. □ 
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2.4- An admissible stress-free-like boundary condition 

The principal motivation to consider the so-called stress-free boundary 
condition is that it minimizes viscous layers and is thus less computationally 
demanding than the no-slip boundary condition. We have seen above that 
this boundary condition unfortunately leads to pathological stability proper- 
ties when the computational domain is axisymmetric. A possible remedy to 
this problem is to consider the following non-symmetric boundary condition: 

(n-Vu) xn|r = 0. (2.9) 

The tangent components of the normal derivative of the velocity field are 
zero. The physical interpretation of this condition is definitely less appealing 
than that of the stress-free boundary condition, but once one realizes that 
the stress-free boundary condition is ad hoc, one comes to think that (2.9) 
is not more ad hoc than the stress-free condition. The main advantage we 
see in (2.9) over the stress-free condition is that it yields standard stability 
properties, i.e., {0} is the global attractor when there is no forcing. 

Lemma 2.3. The following holds for all smooth solenoidal vector field u that 
satisfies (n-Vu) xn| r = 0: 



/ -V-(e(u))-v = \ [ Vu:Vv, Vv G H 1 ^), vn| r = 0. 
Jn Jn 



(2.10) 



Proof. Upon observing that V(e(u)) = |V-(Vu) since u is solenoidal, we 
infer that 

f -V(e(u))-v = f -|V-(Vu).v = \ [ Vu:Vv - \ [ (n-Vu)v 

Jn Jn Jn Jr 

= \\ Vu:Vv-i / (n-Vu)-((n-v)n) + \ / ((n-Vu) xn)- (nxv) 
Jn Jr Jr 



Vu:Vv, 



n 



where we used again the decomposition V|r = (n-v)n — nx (nxv). □ 

Proposition 2.4. Assume that Q is an open, connected, bounded Lipschitz 
domain, then {0} is the global attractor of (2.1)-(2.2)-(2.3)-(2.9)-(2.5). 
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Proof. Repeat the argument in the proof of Proposition 2.1 using Lemma 2.3 
together with the following Poincare-like inequality 



which can be shown to hold by proceeding as in the proof of Proposition 2.2. 



3. Precession driven flow with Poincare stress 

If the fluid domain is a spheroid that undergoes precession, the steady 
state Navier-Stokes equations with the slip condition admit a so-called Poincare 
solution. We show in this section that, independently of the value of the vis- 
cosity, the Poincare solution is not an attractor of the problem if the tangen- 
tial stress at the boundary is enforced to be equal to that of the steady-state 
Poincare solution. 

3.1. Geometry and equations 

The container is an ellipsoid of revolution of center O and symmetry axis 
Oz. The unit vector along the Oz-axis is e z . The unit vectors along the other 
two orthogonal axes Ox and Oy are e x and e y , respectively. The surface of 
the ellipsoid is defined by the equation 



where (5 > —1 and (3 ^ 0. We assume that the container rotates about the 
Oz-axis with angular velocity e z and that this frame slowly precesses about 
the Ox-axis with angular velocity ee x (this particular precession angle is 
investigated in [17]). The non-dimensional Navier-Stokes equations describ- 
ing the motion of the fluid in the non-inertial precessing frame of reference 
(O, e x , e y , e z ) are written as follows: 





□ 



x 2 + y 2 + (l + (3)z 2 = l, 




d t u + u Vu - 2i/V-e(u) + 2ee I xu + Vp 

Vu 



0. 



(3.2) 
(3.3) 
(3.4) 



u . 



We additionally enforce the slip boundary condition, 



u-n|r = 0. 



(3.5) 
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The system (3.2)-(3.3)-(3.5) is known to admit a steady solution called 
the Poincare flow (see e.g. [17]); its expression is: 



Similarly to [17] we consider the problem (3.2)-(3.3)-(3.5) equipped with the 
additional non- homogeneous boundary condition 



That is, we want the tangential component of the normal stress to be equal 
to that of the Poincare solution. As mentioned in [17], it is clear that 

Claim 3.1 (See [17]). u P is a steady state solution of (3.2)-(3.3)-(3.5)-(3.7). 

3.2. Long term stability 

The question that we now want to investigate is whether there is a thresh- 
old on v beyond which u P is a stable solution as t — > +oo; i.e., does the flow 
return to u P independently of the initial data as t — > +oo if v is large 
enough? We show in this section that the answer to this question is no, 
the fundamental reason being that solid rotations cannot be dampened by 
viscous dissipation, no matter how large v is. 

Proposition 3.1. For all v > 0, {up} is not an attractor of the Navier- 
Stokes problem (3.2)-(3.3) equipped with the boundary conditions (3.5)-(3.7). 

Proof. Let p > be an arbitrary positive number. Let B(u P ,p) C H be a 
ball of radius p centered at u P . Let w = ue 2 xr is a solid rotation about the 
Oz-axis, and assume that u ^ is small enough so that u P + w G B(u P , p). 
Let us prove that u P + w is a steady state solution of (3.2)-(3.3)-(3.5)-(3.7). 
Owing to e(w) = 0, w-ri|r = 0, V-w = 0, it is clear that u P + w is solenoidal 
and satisfies the boundary conditions (3.5)-(3.7). Let us now show that it is 
possible to find a pressure field so that the steady state momentum equation 
holds. 

Let us first prove that u P -Vw + w-Vup + 2ee x xw is a gradient. A 
straightforward computation gives: 




(3.6) 



(n-e(u)) xn !r = (n-e(u P )) xri| r . 



(3.7) 
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so that 



up-Vw + w-Vup + 2ee x xw = -V (u(x 2 + y 2 )) + V (^jfi 1 + Z 3 )^ 



Let use then define g(x) := — u(x 2 + y 2 ) + ^(1 + f3)xz. Observe that we 
can define the pressure field r(x) so that Vr := -up-Vup — 2ee x xup, since 
up solves (3.2). Let us finally observe that w-Vw = — |V|w| 2 . Then we 
conclude that up + w solves (3.2) with p = q + r — ^|w| 2 . In particular if we 
set u = up + w, then up + w remains a solution forever, i.e., the solution 
does not converge to Up as t — > +oo, no matter how small p is and no matter 
how large v is. □ 

3.3. Angular momentum balance 

Let us now mention a result on the balance of the angular momentum. 
Let us assume that u solves (3.2)-(3.3) with the boundary conditions 

nu = on T, (3.8) 

(n-e(u)) xn = gxn on T, (3.9) 

where the field g is a boundary data. Let us now define the angular momen- 
tum 

M := / xxu. (3.10) 
Jn 

Lemma 3.1. Denoting by M z and M y the z- and y -component o/M, respec- 
tively, all the weak solutions of (3.2) -(3.3) -(3.8) -(3.9) satisfy 

d t M z + eM y = - [/(gxn)'((e z xx)xn), a.e. t e (0, +oo). (3.11) 
Jda 

Proof. Observing that M z = J^(e 2 xx)-u, we multiply (3.2) by e 2 xx and 
integrate over Q. Using the divergence free condition together with (3.8) and 
integrating by parts, we infer that 

/ (e 2 xx)-(u-Vu)= / V-(u <g) u)-(e z xx) = / (u-n) (u- (e 2 xx)) = 0, 
Jn Jn Jan 

where we used that (u®u):V(e 2 xx) =0 since the matrix u(g)u is symmetric 
and V(e 2 xx) is anti-symmetric. The same argument applies to the viscous 
term 

/ (e 2 xx)-z/V(e(u)) = / z/(e(u)-n)-(e 2 xx) = / i/(gxn)-((e 2 xx)xn), 
Jn Jan Jon 



13 



where we used e 2 xx = (e 2 xx)xn since (e z xx)-n| r = 0. The same argument 
applies again for the pressure term since Vp = V-(p/) where / is the identity 
matrix. 

/ (e 2 xx)-Vp= / p(e 2 xx)-n = 0. 
Jn J 80, 

We now deal with the Coriolis term by applying Lemma 3.2: 
|je 2 xx)-( ej xu) = ^ jT e„-(xxu) = ^M y . 

The conclusion follows readily. □ 

Lemma 3.2. Let v e L 1 (fi) be an integrable vector field such that Vv = 
and v-n| r=0; then 



e s -(xxv) = 2 / (e z xx)'(e x xv). (3.12) 



n 



Proof. Let us first observe that | SJ (e 2 xx)-(e a .xv) = — f Q xu z . Noticing that 
J n xv z + zv x = f n v-V(zx) = since V-v = and v-ri| r = 0, we infer that 

/ (e 2 xx)-(e :E xv) = - / xu z = \ [ zv x -xv z = \ [ e„-(xxv), 
Jn Jn * Jn * Jn 

which concludes the proof. □ 

Remark 3.1. If we choose g = e(up)-n like in (3.7), then — J m z/(gxn)-((e 2 xx) 
is equal to — | S] (e 2 xx)-i/V(e(up)) = and the balance equation of the an- 
gular momentum in the z direction simplifies to d t M z + eM y = 0. 

Remark 3.2. Note that (3.11) is just a consequence of (3.2)-(3.3)-(3.8)-(3.9). 
This balance holds whether the long term stability of (3.2)-(3.3)-(3.8)-(3.9) 
is spurious or not. It is false to consider that (3.11) is an additional equation 
that fixes the long term stability behavior of (3.2)-(3.3)-(3.5)-(3.7). 



4. Precession driven flow with stress-free boundary conditions 

We show in this section that if we enforce e(u)-n| r = 0, instead of en- 
forcing e(u)-n| r = e(u P )-n| r in (3.2)-(3.3)-(3.8), then becomes the unique 
stable solution as t — > +oo, i.e., {0} is the global attractor. 
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4-1. Long time stability 

The setting of the problem is the same as in Section 3.1 except that 
we enforce the tangential component of the normal stress to be zero at the 
boundary. 



2uVe(u) + 2ee x xu + Vp 


= 


in Q 


(4.1) 


Vu 


= 


in Q 


(4.2) 


nu 


= 


on r 


(4.3) 


(n-e(u)) xn 


= 


on r 


(4.4) 


U|t=o 


= u 


in Q. 


(4.5) 



The result that we want to emphasize is that contrary to what we observed 
in Section 3, becomes the unique stable solution of (4.1)-(4.5) as t — > +oo. 
The main result that we want to prove here is that any solution of the 
system (4.1)-(4.4) returns to rest as t — > +oo. The key argument is that 
solid rotations about the Oz axis are not stationary solutions of (4.1). This 
fact has been mentioned in [17] without proof. 

Theorem 4.1. {0} is the global attractor of (4.1)-(4.5). 

Proof. Let us start by observing that {0} is indeed an invariant set of (4.1)- 
(4.5). Let B(0,p) be the unit ball in H centered at and of radius p. Let 
u G B(0,p) and let u G L 2 ((0, +oo); L 2 (fi)) n L°°((0, +oo); H 1 ^)) be a 
Leray-Hopf solution of (4.1)-(4.5) and consider the following decomposition: 

u(t) = ir'-(t) + A(t)e 2 xx, where u x (t) G TZ ± , X(t) el, Vt G [0, +oo). 

Lemma 2.2 together with u being a Leray-Hopf solution implies that 

W^mh^+lHtf + ^K [ llu^H^njdr < \\u \\l Hn y 

Jo 

where 7 = ||e 2; xx|| 2 j2 . Using the Gronwall-Bellmann inequality, we infer that 

11^(011^(0) < \\^ohnn)e- 2uKt . 

Let t 2 > ti in (0, +00), then (3.11) means that 

(A(* 2 ) - A(*i)) 7 = -e I" [ e J/ -(xx(A(r)(e,xx) + u ± )). 

Ju Jo, 
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But Lemma 3.2 implying that 

/ e y -(xx(A(r)(e z xx))) = A(r) / e 9 -(xx(e 2 xx)) 
Jn Jn 

= 2A(r) / (e z xx)-(e x x(e z xx)) = 0, 
Jn 

we finally infer that 

|A(i 2 ) - A(*i)| < T X e f 2 [ le^xxu^l < c (e~ 2 ^ - e" 2 ^ 2 ), 

Ju Jn 



where c is a generic constant that depends on Q, u, and p and may vary at 
each occurrence from now on. Note in passing that this also proves that X(t) 
converges to a real number Aoo as t — )■ +00, and |Aoo — A(t)| < ce~ 2uKt . 

Let us take cp G T>(Vi) independent of time and divergence- free. Since 
u is a Leray-Hopf solution (recall t 1 — > u(t) is continuous in the L 2 -weak 
topology) we have 

"*2 



0= / (u(t 2 ,x) -u(t 1 ,x))-^(x)dx 



2eu(t, x)-(<^(x)xe :c ) dx dr 



*2 



ti Jn 



2z/u(r, x)-V-e(y>) dxdr — 



h Jn 



ti Jn 



'u(x) <g> u(x)):V<^(x)dxdr. 



Let us now set t 2 = + Upon observing that J Q ((e 2 xx)®(e 2 xx)):V<^(x)dx 
and J n (e z xx.)-V-e(cp) dx = 0. This implies that there is a constant c(cp) > 
so that 

A(r)dr f (e 2 xx)-(^(x)x ei )dx 
Jn 

+ c(^)e- 2 ^. 

Let us choose (p so that 2e J SJ (e z xx)-(^(x)xe ;i .) dx = 1. The above estimate 
implies that 



2s 



< 



(A(t 2 ) - A(fi)) / (e 2 xx)- ¥ >(x)dx 



/■ti+i 
/ 

Jti 



A(r) dr 



< c(y)e" 



-2vKt 



This in turn implies that A^ = lim^oo J" 1 A(r) dr = 0, which means 
Ano = 0. In conclusion 



which concludes the proof. 



lim ||u(f)|| L 2 ( m < c(<p) lim e" 



-2vKt 



0, 



(4.6) 

□ 
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5. Numerical illustrations 

To illustrate the above mathematical results, we have performed two series 
of numerical simulations similar to those presented in [17]. The authors study 
therein the dynamo action in an oblate spheroid defined by equation (3.1) 
with (5 = 0.5625 (this corresponds to the value b = 0.8 for the semi-minor axis 
used in [17], b := (1 + f3)~^). This spheroid rotates about the Oz-axis and 
precesses about the Ox-axis. Two sets of boundary conditions are considered: 
either the homogeneous stress-free boundary or the Poincare stress condition 
is enforced. Simulations are carried out using a mixed Fourier decomposition 
and finite element code described in details in [6]. 

The first simulation solves the equations (4.1)-(4.2)-(4.3)-(4.4) with the 
initial data u\ t=0 = 0.1(—ye x + xe y ). The precession rate is e = 0.25 and the 
reciprocal of the viscosity is 1/V = 0.024. Figure 1 shows the time derivative 
of the total energy Ek = |||u|Il2 m the precessing frame. Note that d t Ex is 
always negative, establishing that Ek is a decreasing function. This graph 
is in excellent agreement with figure 1 of [17]. It also shows that u — > as 
t — > oo in agreement with (4.6) (i.e., {0} is indeed the attractor) . 




" 5 20 40 60 80 100 120 
time 



Figure 1: Time evolution of d t EK of the solution of equations (4.1)-(4.2)-(4.3)-(4.4) for 
/3 = 0.5625, e = 0.25 and l/v = 0.024. 

The second series of simulations solves equations (3.2)-(3.3)-(3.5)-(3.7) 
with an ad hoc initial condition which is a very small perturbation of the 
Poincare solution. The parameters are e = 0.25 and l/v = 0.00375. Figure 2 
shows the time evolution of the kinetic energy of the perturbation to the 
Poincare solution, 5E^ = |||u — Up||^2^), from t — to t — 1100, see curve 
labeled "0 perturb" . The energy grows exponentially initially, then saturates 
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around an oscillatory state. These results are similar to those shown in figure 
1 of [17]. 

In order to evaluate the influence of solid rotations, we restart the com- 
putation at t — 1100 by adding the perturbation ±0.025(— ye x + xe y ) to 
the solution. Since the maximum norm of the Poincare solution is 1.25, the 
added perturbations are only 2% of the maximum velocity. Time integration 
is performed in each case until convergence to an oscillating state is obtained. 
The curves corresponding to the time evolution of the kinetic energy of the 
solutions thus obtained are labeled "0.025 perturb" and "—0.025 perturb" 
in Figures 2-2 (b). We observe in Figure 2(b) that these perturbations have 
strong impacts on the asymptotic solutions. To better compare our results 




600 700 800 900 1000 1100 1200 ' 1050 1100 1150 1200 1250 



time time 

(a) 5E K vs. time (b) zoom 

Figure 2: (Color online) Time evolution of the kinetic energy, SEk, of the perturbation 
of the solution of (3.2)-(3.3)-(3.5)-(3.7) with /3 = 0.5625, e = 0.25 and 1/v = 0.00375 (a), 
and zoom (b). 

with those from [17], we show in figure 3 the energies in the northern 5Exn 
and southern 5E Ks hemispheres of the spheroid. The perturbation with the 
positive sign increases both the total kinetic energy and the amplitude of the 
northern and southern energies, whereas the perturbation with the negative 
sign decreases both the total kinetic energy and the amplitude of the oscilla- 
tions of the northern and southern energies. The oscillations of the northern 
and southern energies obtained with the positive perturbation are more sinu- 
soidal than those obtained with the negative perturbation. The shape of the 
oscillations of the northern and southern energies obtained with the negative 
perturbation are similar to those in [17] (see the more pronounced nonlinear 
shape). These simulations illustrate well that the 0-perturbation solution is 
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Total - 
North - 
South .v 



(a) SE K , 5E Kn , 8E Ks with 
ad hoc initial condition 



(b) 0.025 perturb 



(c) -0.025 perturb 



Figure 3: (Color online) Kinetic energy, SEk = |||u — upH^m, where u solves (3.2)- 
(3.3)-(3.5)-(3.7) with f3 = 0.5625, e = 0.25 and 1/v = 0.00375: on each graph, top curve 
is SEk, bottom curves are the energies SEx n (dashed line) and 8Ek s (dotted line) in the 
northern and southern hemispheres. 



not an attractor, i.e., it is not stable under perturbations. We have verified 
(results not shown here) that an entire family of solutions can be obtained 
from the 0.025-perturbation to the — 0.025-perturbation solutions by scaling 
the perturbation appropriately. These tests show that using the stress-free 
boundary condition to evaluate nonlinear behaviors of Navier-Stokes systems 
may sometimes be dubious when the domain is axisymmetric. 



6. Discussion 

The so-called stress-free boundary condition (n-e(u))xn| r = is often 
used in the geodynamo literature to avoid issues induced by viscous layers. 
For example, very recently, an anelastic dynamo benchmark [7] was con- 
ducted in a rotating spherical shell. The authors emphasize in their conclud- 
ing section the difficulties they encountered to compare four different codes 
using a model with stress-free boundary conditions applied to the ICB and 
the CMB. Since the container is a spherical shell, the balance equation (3.11) 
gives <9tM = 0, and each group had to apply some remedy in order to numer- 
ically conserve the three components of the angular momentum. But, more 
importantly, they also had to use the same initial condition. There was no 
such difficulties in the older dynamo benchmark [2] using the same geome- 
try because the no-slip boundary condition was prescribed at each interface. 
These results illustrate again that the stress-free boundary condition induces 
pathological stability behaviors when the flow domain is axisymmetric. 

We have shown in this work that stress-free boundary condition leads to 
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spurious behaviors when the fluid domain is axisymmetric. We hope that 
the present work will help draw the attention of the geodynamo community 
on this problem. The above pathological stability behaviors can be avoided 
by enforcing one additional condition. For instance, for problem (3.2)-(3.5) 
and (3.7), one could think of enforcing the vertical component of the angular 
momentum of the perturbation to the Poincare flow, say 

(u — up)-(e 2 xx) ds = 0, (6.1) 

or enforcing the perturbation and the Poincare flow to be orthogonal in 
average over the boundary, say 

(u-up)-upds = 0. (6.2) 

For problem (2.1)-(2.5), one could think of enforcing the vertical component 
of the total angular momentum 

u-(e 2 xx) ds = 0, (6.3) 

as was done for the three components in the anelastic dynamo benchmark [7] . 

We have suggested in §2.4 to use a boundary condition that does not 
have the stability problems mentioned above. For the problem (3.2)-(3.5) 
this condition is 

(n-Vu)xn| r = (n-Vu P )xn| r , (6.4) 
and for the problem (2.1)-(2.5) this condition is 

(n-Vu)xn|r = 0. (6.5) 

Let us finally emphasize that it is false to consider that the momentum 
balance equation (3.11) is an additional equation that makes (3.2)-(3.3)-(3.5)- 
(3.7) a well-behaved dynamical system. The equation (3.11) is a redundant 
consequence of (3.2)-(3.3)-(3.5)-(3.7). For instance, (6.1) (or (6.2) or (6.3)) 
is an additional equation whereas (3.11) is not. 
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